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ABSTRACT 

Magnetic perturbations to the frequencies of low degree, high radial order, ax- 
isymmetric pulsations in stellar models permeated by large scale magnetic fields are 
presented. Magnetic fields with dipolar, quadrupolar and a superposition of aligned 
dipolar and quadrupolar components are considered. The results confirm that the 
magnetic field may produce strong anomalies in the power spectra of roAp stars. It is 
shown for the first time that anomalies may result both from a sudden decrease or a 
sudden increase of a mode frequency. Moreover, the results indicate that the anomalies 
depend essentially on the geometry of the problem, i. e. on the configuration of the 
magnetic field and on the degree of the modes. This dependence opens the possibility 
of using these anomalies as a tool to learn about the magnetic field configuration in the 
magnetic boundary layer of pulsating stars permeated by large scale magnetic fields. 
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1 INTRODUCTION 

Rapid oscillations, with amplitudes of a few mmag have 
first been discovered in cool, chem ically peculiar magnetic 
Ap stars over two decades ago bv iKurta <ll982t) . The fre- 
quencies, v, of the observed oscillations are typically within 
the range lmHz < v < 3mHz, a lthough lower frequency 
pulsations hav e been predicted b y lCunhal <l2002t) - and re- 
cently found bv lElkin et. all (120051) - among the more evolved 
cool Ap stars. Recent observational and theoretical reviews 
th is class of pulsator s . kno wn a s roAp stars, are p rovided 
inlKurtz fc Martinez! feOOCft an d lkurtzl \2004) . and ICunhal 
J2005ll2005l) and lGouehl i2005l) . 

The oscillations observed in roAp stars are believed to 
be high order acoustic modes modified near the surface by 
the magnetic field. In principle, the high order of the ob- 
served oscillations should simplify the theoretical analysis 
of the pulsation power spectra of roAp stars, by allowing 
the use of standard asymptotic tools. For this reason, roAp 
stars have been recognized as excellent potential targets for 
asteroseismology . How ever, as emphasized in t he works of 
iMatthews et all (II999T) and lCunha et alJ J2003lh the pulsa- 
tion power spectra of these stars cannot be reconciled with 
that produced with standard stellar and pulsating models. 
Part of the discrepancy may have its origin in the difficulty 
of deriving accurate effective temperatures and bolometric 
corrections (a nd hence radii) for Ap stars (lLeckronelll973t 
IStepierilll99lh However, the complexity of the surface lay- 
ers of these stars, including the presence of magnetic fields 



and chemical inhomogeneities indicates that a more sophis- 
ticated modeling of the pulsations is also desirable when 
attempting to interpret their power spectra. The present 
work builds on previous studies of the direct effect of the 
magnetic field on pulsations in the presence of large-scale 
magnetic fields. A detailed discussion of the impact that the 
results presented here may have on the power spectra of 
selected models of roAp stars will be provided in a future 
paper. 

1.1 Direct and indirect effects of the magnetic 
field 

1 The strongest magnetic field so far detected in a cool 
Ap star has a mean magn etic field modulus of < B > = 
25.5kG jHubrig et al.ll2005l) . Generally speaking, however, 
roAp stars are permeated by only relatively strong, large 
scale magnetic fields, with typical magnitude s of a few kG 
jMathvs fc Hubrid Il997l: iHubrig et alJ 120041) . These mag- 
netic fields influence the oscillations both directly, by gen- 
crating an additional restoring force that affects the wave 
dynamics in the surface layers, and indirectly, by interfering, 
and possibly suppressing, envelope convection. In fact, the 
suppression of envelope convection, at least in some angular 
region of the star, seems to be a necessary condition for the 
driving of the observed high frequency osci llations in other- 
wise standard stellar models o f roAp stars fealmforth et all 
l200lt ICunhall2002l: ISaidl2005D . 

Magnetic perturbations to the eigenfrequencies and 
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eigenfunctions of the oscillations in stellar models perme- 
ated by an intense, large scale magnetic field have been com - 
puted bvlDziembowski fc Goodd (llOQfiTl.lBigot et al l l)2000Fl . 
ICunha fc GougblfcOOOl) and lSaio fc Gautschvl (120041) . These 
authors used a singular perturbation approach, that takes 
into account the fact that in the surface layers the magnetic 
effect on t he oscillations cann ot be treated as a small per- 
turbation feiront et alJll982l) . These works revealed some 
important consequences of the direct effect of the magnetic 
field on pulsations. In particular, it was found that the eigen- 
functions are perturbed to an extent that at the surface they 
can no longer be described by a single spherical harmonic 
function. When the cancellation effects originated by our 
incapacity to resolve the star are taken into account, it is 
clear that this can lead to an erroneous identification of the 
degree of the mode that is observed. In fact, if excited, it is 
possible that modes of degree significantly higher than I = 3 
may generate lower degree components near the surface that 
may be observed. 

Another important result was the finding that for cer- 
tain magnetic field intensities and oscillation frequencies the 
pulsation energy losses are particularly high. These energy 
losses may contribute to the selection effects that determine 
the frequency range in which unstable high frequency oscil- 
lations may be found, for a given magnetic field. Finally, the 
maxima of energy losses mentioned above were also found 
to be associated with abrupt changes in the oscillation fre- 
quencies. These abrupt changes may introduce anomalies in 
the power spectra of roAp stars, as was suggested to hap- 
pen i n one of the best studied roAp stars HR1217 dCunhal 
2001). Such abrupt changes will be discussed in more detail 
in sect ion El 

In lCunha fc Gougbl i2000T) a variational approach to the 
calculation of the frequency shifts was applied to polytropic 
models of roAp stars. In the present paper we extend their 
work, applying the same approach to standard stellar mod- 
els permeated by magnetic fields of different configurations. 
As before, only the direct effect of the magnetic field on pul- 
sations is taken into account. In section |2] we set the prob- 
lem and discuss the physics of the pulsations under study 
in the different regions of the star. In section |H] we describe 
the equilibrium model and the magnetic field configurations 
that will be considered. 

In section 0] we recall the expression derived for the 
magnetic freque ncy perturbations under the variational ap- 
proach used by ICunha fc Gougbl i2000F) and in section 
we discuss the results obtained in particular for a) different 
boundary conditions applied to the pulsations near the sur- 
face and b) different magnetic field configurations. A com- 
parison with the results obtained by Saio (private communi- 
cation) for the same model is also presented in that section. 
Finally, in section [5] we summarize our most relevant results. 



2 MAGNETOACOUSTIC WAVES 
2.1 Solutions 

Adiabatic motions in a non-rotating plasma in the limit of 
perfect conductivity are described by the system of magne- 
tohydrodynamic equations, 
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and no is the permeability of the vacuum, p is the density, 
p is the pressure, g is the gravitational field, B is the vec- 
tor magnetic field, v = |j| is the velocity, £ is the vector 
displacement, 7 is the first adiabatic exponent and t is the 
time. 

In a star permeated by a magnetic field, the direct con- 
tribution of the field to the dynamics of what would other- 
wise be acoustic pulsations is important only near the sur- 
face, where the magnetic pressure is comparable, or larger, 
than the gas pressure. We shall refer to that region as mag- 
netic boundary layer. 

The system of equations Q-JSJ is solved, for the mag- 
netic boundary layer, in a local, plane-parallel approxima- 
tion, assuming small perturbations to the equilibrium struc- 
ture, which is assumed to be permeated by a force free 
magnetic field. The (Eulerian) perturbation to the gravi- 
tationl potential is ignored (Cowling approximation). More- 
over, only high order, low degree, axisymmetric modes will 
be considered. The equations that are actually solved are 
described in appendix A. In the present section we attempt 
only to give some insight into particular classes of solutions 
that might be found under the conditions that are of interest 
to us. 

Let us consider first the simple case of a homogeneous 
fluid, permeated by a uniform magnetic field, in conditions 
such that the magnetic and gas pressure are comparable. 
Under these conditions, two types of linear, adiabatic, non- 
diffusive magnetoacoustic perturbations may be identified, 
namely, the fast and the slow magnetoacoustic waves (e.g. 
|Priestm982Tl . These waves correspond to two distinct solu- 
tions of the equations that govern the perturbations in such 
a homogeneous fluid. They may be regarded as an acoustic 
wave modified by the magnetic field and as a compressional 
Alfven wave modified by the gas pressure. 

When the density is stratified, as in a real star, it is gen- 
erally no longer possible to identify the two distinct solutions 
described above. Instead, the perturbed equations describe 
a single, rather more complex, magnetoacoustic wave. How- 
ever, in the outermost layers of the star, where the gas pres- 
sure is very small compared to the magnetic pressure, and 
below the surface layers, where the gas pressure dominates 
over the magnetic pressure, the magnetoacoustic wave de- 
couples into waves that are essentially magnetic and waves 
that are essentially acoustic. Thus, while the equations de- 
scribing the linear, non-diffusive, magnetoacoustic perturba- 
tions have to be solved numerically throughout the region 
where the magnetic and gas pressure are comparable, ap- 
proximate analytic solutions can be found to describe the 
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Table 1. global parameters of the model used in the computa- 
tions. M/Mq and _R/Rq are the mass and radius in solar units, 
T e ff is the effective temperature and log(L/LQ) is the logarithm 
to the base 10 of the luminosity in solar units. 



Model 


M/M 


R/Rq 




log(L/L Q ) 


Cesam - ZAMS 


1.8 


1.57 


8362 K 


1.035 



decoupled acoustic and magnetic waves both in the upper 
atmosphere and in the interior of the star. 

2.2 Boundary conditions 

Some of the boundary conditions applied in the present work 
to solve the perturbed equations in the magnetic boundary 
layer make use of the decoupling described above. In partic- 
ular the numerical solutions are matched into approximate 
analytical solutions both in the outer layers of the star and 
in the interior. The details of these and of the remaining 
boundary conditions applied are given in appendix B. Never- 
theless, it is worth noting here that two possibilities are con- 
sidered regarding the reflection/transmission of the waves 
in the surface layers. In some cases a mechanical boundary 
condition is applied which assures that the waves are fully 
reflected at the outermost layer of the model, regardless of 
the oscillation frequency. In other cases, the component of 
the displacement parallel to the magnetic field direction in 
the atmosphere is matched onto an outgoing propagating 
wave whenever the oscillation frequency is above the acous- 
tic cutoff frequency. This is justified by the fact that in the 
upper atmosphere the component of the solution parallel to 
the magnetic field direction is essentially an acoustic wave 
(Sousa & Cunha - private communication) and, thus, its en- 
ergy is expected to be lost when the oscillation frequency 
is above the acoustic cutoff frequency. We recall that the 
acoustic cutoff frequency appropriate for the oscillations un- 
der study depends on the inclination of the magnetic field, 
being equal to the cut-off frequency in the absence of the 
magnetic field when the direction of the latter is locally ver- 
tical and tending to zero when th e direction of the local 
magnetic field tends to horizontal jDziembowski fc Goodd 
ll99fj:lBigot et al.ll200d) . 



3 MODELS 

The equilibrium model used in this work was compu ted by 
Marq ues et al. 1 with the evolutionary code CESAM iMorell 
1997) . The atmosphere of the model was extended to a 
minimum density of 10 _11 <?cm~ 3 . The global parameters of 
the model are given in tabled 

Since the magnetic field is assumed to be force free, 
it does not influence the stratification of the equilibrium 
model. However, when small perturbations to the equilib- 
rium structure are considered, the perturbed magnetic field 

1 models available at www.astro.up.pt/corot/models/cesam/ 



responds through the Lorentz force, j A B, present in equa- 
tion ©. Thus, we shall refer to 'the magnetic model' and 
'the non-magnetic model' whenever we compare the oscil- 
lations in the presence of a magnetic field with the oscil- 
lations in an otherwise similar model without a magnetic 
field. Moreover, the terms perturbed and unperturbed will 
hereafter be used when referring to the properties of the 
oscillations in the magnetic and non-magnetic models, re- 
spectively. 

The influence that the magnetic response has on what 
would otherwise be an acoustic oscillation, depends, among 
other factors, on the configuration of the magnetic field. 
Three different magnetic field configurations are considered 
in the present work: a dipolar and a quadrupolar magnetic 
field defined, respectively, by 

B d = 6 d /r 3 (cos6» e r + 1/2 sin 6 eg), 

and 

B q = 6 q /r 4 (l/2(3cos 2 6i - 1) e r + cos 6 sin 6 e e ), 

and a magnetic field, Bd q , composed of aligned dipolar plus 
quadrupolar components. Here fed and 6 q are constant am- 
plitudes, r and 8 are the usual coordinates in the spherical 
coordinate system (r, 9, <f>) and e r and ee are the correspond- 
ing unit vectors. 

Since the magnetic perturbations involve only quadratic 
terms in the magnetic field, the perturbations induced on the 
eigenfunctions by dipolar or quadrupolar magnetic fields are 
necessarily symmetric about the magnetic equator. Thus, 
when expanded in spherical harmonics, the perturbed eigen- 
functions will include only degrees I of parity identical do the 
degree of the unperturbed mode. On the other hand, the ef- 
fect of a magnetic field whose configuration is the sum of 
dipolar and quadrupolar components does not comply with 
the referred symmetry properties. Consequently, unlike the 
dipolar and quadrupolar cases, the third magnetic field con- 
sidered here can perturb the displacement eigenfunctions in 
a way such that their expansion in spherical harmonics con- 
tains degree components of different parity. 



4 MAGNETIC PERTURBATIONS: 
VARIATIONAL APPROACH 

The direct effect produced by the magnetic field on the wave 
dynamics in the outer layers of the star modifies the frequen- 
cies of the oscillations. We define the magnetic frequency 
perturbations as the difference between the eigenfrequencies 
appropriate to the magnetic and the non-magnetic models. 
Different approaches to the calculation of these frequency 
perturbations ar e possible. Here we foll ow the variational 
approach used in lCunha fc Gouehl i2000t) . 

In the interior of the magnetic model the magnetic field 
has a negligible contribution to the restoring force associated 
with the fast wave. There, the dynamics is effectively field- 
free and the acoustic oscillations are described by the usual 
equations of hydrodynamics. Consequently, when compared 
with the acoustic wave propagating in the interior of the 
non-magnetic model, the fast wave in the interior of the 
magnetic model differs only due to the effects introduced by 
the overlying magnetic boundary layer. 

In the asymptotic limit valid for modes of high radial 
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order, the radial component of the displacement in the non- 
magn etic model i s given, in the region of propagation, by 
fe.g. lGouglJll99fl 



(r, 9, 4>,t) ~ A 



k 1 ' 2 

1/2 

rPo 



ndr + 5 I Yj m e 



(7) 



where po is the density in the equilibrium model, k is the 
vertical acoustic wavenumber, u is the oscillation frequency, 
A is a constant, 8 is a phase, R* is a particular value of the 
radial coordinate and Y ; m is a spherical harmonic of degree 
I and order m. 

In the magnetic model, the spherical symmetry of the 
problem is broken by the effect produced by the magnetic 
field on the oscillations in the magnetic boundary layer. Con- 
sequently, the dependence on latitude of the radial displace- 
ment in the interior of the magnetic model is no longer given 
by a single spherical harmonic of degree I. However, the fact 
that the magnetic field varies only on large scales and that 
the modes to be considered are of low degree allows the 
problem to be solved locally, assuming at each latitude a lo- 
cally uniform field (see appendix A for details). Under this 
assumption the radial component of the displacement be- 
low the magnetic boundary layer (but still sufficiently close 
to the surface), may still be expressed, at each particular 
latitude, by equation J7J. 

The effect of the magnetic boundary layer on the oscil- 
lations in the interior may be seen as a shift in the phase S 
computed at the base of that layer (defined to be r — R*), 
which changes from the value 5 unp it would have in the non- 
magnetic model to a new value 5 which depends on latitude. 
Moreover, the vertical acoustic wavenumber and the oscilla- 
tion frequency are also modified. 

In the variational approach the phase shifts AS (9) = 
5 — <5 un p are used to determine the magnetic perturbations 
to the eigenfrequencies. In practice the variational princi- 
ple is applied to determine the perturbations to the eigen- 
frequencies without having to dete rmine the perturbations 
to the eigenfunctions. According to lCunha fc Gougbl (|2000F) 
the difference between the angular frequencies in the mag- 
netic and non-magnetic models, for a mode of degree I and 
azimuthal order m, is given by: 



Aw 



AS 



f* co -2 K -1 dr' 



(8) 



where Co is the sound speed in the equilibrium model, n 
is the lower turning point (at which « = 0) and AS is the 
integral phase shift, i.e. 

(9) 



Jo(Yi m ) 2 sin6>d0 

From this expression we can calculate the shift in the cor- 
responding cyclic frequency through the relation Av — 
(27r) _1 Aw. The base of the magnetic boundary layer, R*, 
(at which the solutions appropriate to the magnetic bound- 
ary layer are matched onto those appropriate to the interior) 
must be placed in the region where the gas pressure dom- 
inates over the magnetic pressure, but sufficiently close to 
the surface for the plane-parallel approximation to be valid. 
In the present calculations we have used R* / R = 0.98 

At certain frequencies the eigenfunctions might be sig- 



nificantly perturbed by the magnetic field, at least in some 
angular region of the star (cf. works mentioned in subsection 
II. H . This fact raises some concern regarding the accuracy of 
the results obtained for those frequencies, independently of 
which approach (variational or expansion in spherical har- 
monics) is used. 

If the eigenfunctions are strongly perturbed in a rela- 
tively small range of latitu des, the spherical harmonic ex - 
pansion use d in th e worksof jDztembowski fc' Goodd <199fl . 
iBieot et alJ (l2000l) and ISaio fc Gautschvl feOojT mav be 
'contaminated' by rather high I terms. These terms might 
be completely missed if a simple conve rgence criteria is 
used. Such problems were acknowledge in lSaio fc Gautschvl 
(2004), who reported convergence problems at some values 
of the unperturbed frequency. 

In the case of the variational approach, the unper- 
turbed eigenfunctions in the interior are used to calcu- 
lat e the perturbations to the eigenfrequencies. As argued 
bv lCunha fc Gougbl i200Ct) . even though the eigenfunctions 
might be strongly modified in the outer layers, in the in- 
terior, where most of the inertia of the modes resides, the 
perturbation to the eigenfunctions is likely to be sufficiently 
small to justify the use of a first order variational approach. 
However, the possibility that the sharp variations found by 
ICunha fc Goughl i200(Jl in the phase shifts might trap the 
modes in particular angular regions of the star (Montgomery 
& Gough - private communication) raises some worries. If 
this trapping takes place the eigenfunctions in the interior of 
the magnetic and non-magnetic models will be significantly 
different. With this in mind Montgomery & Gough are in- 
vestigating on the possibility of determining the magneti- 
cally perturbed eigenfunctions without recourse to a spher- 
ical harmonic expansion. In a future work we plan to use 
their results to test the accuracy of the frequency shifts de- 
rived with the variational approach applied here as well as 
to determine the corresponding higher order corrections to 
the latter. 



5 RESULTS AND ANALYSIS 

5.1 Effect of boundary condition at the top 

FigurcsQand|2]compare, respectively, the frequency and the 
phase shifts that are obtained with a dipolar magnetic field, 
when two different boundary conditions are applied to the 
acoustic component of the wave at the outermost layer of 
the star. For the phase shifts only the range of co-latitudes 
< 6 < 7r/2 is shown since for magnetic fields with a dipo- 
lar configuration the perturbations are symmetric about the 
magnetic equator. 

In figurc^it is seen that the boundary condition applied 
might be particularly important for the determination of the 
imaginary part of the oscillation frequency, if the frequency 
of the mode is above the acoustic cutoff frequency appro- 
priate to the non-magnetic model, which in the present case 
is v c » 2550/iHz. Comparing the results for two different 
magnetic field intensities (left and right hand side panels of 
figure we also see that the boundary condition applied 
has a larger impact on the results when the imaginary part 
of the frequency shifts is close to a maximum than when it 
is close to zero. 
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Figure 1. Comparison of the real (top) and imaginary (bottom) parts of the frequency shifts, as function of the cyclic frequency v, 
when two different boundary conditions at the outermost layer are applied: full reflection (filled symbols) and partial transmission (open 
symbols) at the top. Results are shown for two magnetic field polar intensities, namely B p = = lkG (left) and B p = = 2kG (right). 
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Figure 2. Comparison of the real (top) and imaginary (bottom) parts of the phase shifts when two different boundary conditions at the 
outermost layer are applied for two different oscillation frequencies. The continuous line shows the case of full reflection and the dashed 
line shows the case of partial transmission at the top. Results are shown for a dipolar magnetic field with polar intensity B p = b d = 1KG, 
for two different unperturbed frequencies, namely, v = f789 /tHz (left) and v = 3094 /tHz (right). The cutoff frequency for acoustic 
modes in the absence of a magnetic field is v c f» 2550 /iHz for the model used. 
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For the phase shifts, two oscillation frequencies are 
shown, namely, below (left panels) and above (right panels) 
v c . When the oscillation frequency is below v c the phase 
shifts are almost independent of the boundary condition ap- 
plied. This was to be expected near the magnetic poles. 
There, the magnetic field is nearly vertical, and the criti- 
cal frequency for the acoustic component of the wave (that 
may be identified in the outer magnetically dominated layer) 
is similar to v c . Thus, for frequencies below v c this compo- 
nent of the wave is reflected independently of which bound- 
ary condition is applied and, as a result, the phase shifts 
obtained are very similar. Away from the magnetic poles, 
however, the critical freque ncy for the acousti c component 
of the wave decreases (e.g. iBigot et al.ll200oT) . tending to 
zero as the magnetic field direction becomes closer to hori- 
zontal. The fact that the phase shifts are very similar also for 
higher co-latitudes indicates that the acoustic component of 
the wave propagating away in the outer layers has a rela- 
tively small effect on the acoustic oscillation in the interior 
when the magnetic field is significantly inclined. 

When the oscillation frequency is above v c a significant 
difference is found between the phase shifts obtained with 
the two different boundary conditions. This difference is sig- 
nificant up to co-latitudes of about 60°, becoming negligible 
as the magnetic field direction becomes closer to horizontal. 

5.2 Different magnetic field configurations 

Figures |3] and 0| show, respectively, the frequency shifts for 
different magnetic field configurations, all with the same po- 
lar strength, and the frequency shifts for modes of differ- 
ent degree. The frequency shifts follow a similar trend in 
all cases: as the radial order increases, the frequency per- 
turbations generally increase smoothly. At some particular 
frequencies, however, the real part of the frequency pertur- 
bation jumps by an amount that is significantly larger than 
the mode-to-mode variations seen in the smooth trend. At 
the same unperturbed frequencies, the imaginary part of the 
perturbation reaches a maximum. This behavior is similar t o 
that first found in this context bv lCunha fc Goughl |2000), 
for polytropic models, a nd later confirmed for m ore sophis- 
ticated stellar models bv lSaio fc Gautschvl i2004T) . 

There are some additional facts that become evident in 
the present results, namely, 

• the jumps in the frequency shifts may be positive or 
negative, i.e., they may increase or decrease the frequency 
of the modes in comparison to what would be expected if 
the smooth trend were to be followed; 

• the frequencies at which the first maxima of energy 
losses take place do not depend on the degree of the modes 
nor on the magnetic field configuration (at least for the mag- 
netic field configurations considered here); 

• the amount by which the real part of the frequency 
shifts change at the jumps depends both on the magnetic 
field configuration and on the degree of the modes. 

Figures |S] and |5] show the real part of the phase shifts 
obtained for two modes of consecutive radial order in mod- 
els permeated, respectively, by a quadrupolar magnetic 
field and a magnetic field composed of aligned dipolar and 
quadrupolar components. 

Each figure has two panels and in each of the panels the 



two modes shown are chosen such that their frequencies are, 
respectively, below and above a frequency jump. It is readily 
seen from these figure that the sudden changes in the real 
part of the frequency perturbations result from jumps of 
~ 7r in the real part of t he phase shifts, which oc cur at some 
well defined co-latitudes. ICunha fc Goughl i200CJ) noted that 
for polytropic models with a dipolar magnetic field the co- 
latitude at which the real part of the phase shifts jumped was 
approximately the co-latitude at which the product B X B Z 
was maximal. For magnetic fields with dipolar configuration 
that co- latitude is 9 — 45° and corresponds to the co- latitude 
at which \B X \ = 1/2|B«|. 

Figures |5] and HJ show that the latter relation holds also 
for magnetic fields with different configurations, in more so- 
phisticated stellar models. Because a quadrupolar magnetic 
field satisfies the condition \B X \ = 1/2\B Z \ twice in the range 
of co-latitudes < 9 < 7r/2, two different jumps in the real 
part of the phase shifts are found, one at 6 « 24° and one 
at 8 ~ 78°, at different oscillation frequencies. These phase 
jumps are responsible for the two frequency jumps that may 
be seen in the frequency range 1000/iHz < v < 2000/iHz in 
the upper right hand side (rhs) panel of figure 01 The fact 
that the second frequency jump is apparent only for modes 
of even degrees ( I — (circles) and I = 2 (squares)) is a 
direct consequence of the high value of 8 at which the real 
part of the phase shift jumps in this case. Modes of odd de- 
gree hardly feel the second jump in the phase shift because 
their amplitudes are very low at the co-latitudes at which 
the phase shifts of the two consecutive modes is significantly 
different. 

In the case of a magnetic field composed of aligned dipo- 
lar plus quadrupolar components the phase shifts are not 
symmetric about the magnetic equator. Hence, unlike the 
case of dipolar and quadrupolar magnetic fields, for this con- 
figuration jumps in the phase shifts may take place in each 
hemisphere separately. In figure |S| it is possible to see the 
first two jumps in the phase shifts for the case of a magnetic 
field composed of aligned dipolar and quadrupolar compo- 
nents of the same polar strength. The first jump in the phase 
shift takes place at 9 w 31.9° while the second takes place 
at 6 » 100.3°. In both cases the condition \B X \ — 1/2\B Z \ 
is satisfied. These jumps in the phase shifts are responsi- 
ble for the two frequency jumps seen in the squared sym- 
bols of the upper panels of figure |3J in the frequency range 
IOOOajHz < v < 2000/iHz. 



5.3 Cyclic behaviour and scaling 

ICampbell fc Papaloizou| dl986l) noted that the equations de- 
scribing the wave dynamics in the magnetic boundary layer 
of a polytropic model are invariant under a particular scal- 
ing of the variables involved. As argued in lCunha fc Goughl 
( 2000), this scaling implies that in polytropic models the fre- 
quency shifts obey the relation Au ~ f{Bp^v), where / is 
some function and B p is the polar strength of the magnetic 
field. In practice, this means that the frequency perturba- 
tions scale in a very particular way with the magnetic field 
and the oscillation frequency . 

ICunha fc Goughl (l2000f) also found that in polytropic 
models the magnetic perturbations vary cyclically with the 
oscillation frequency. Given the scaling mentioned above, 
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Figure 3. Real (top) and imaginary (bottom) parts of the frequency shifts, as function of unperturbed frequency, for magnetic fields 
with polar strength of B p = 3 kG and different configurations - dipole (triangles), quadrupole (diamonds), aligned dipole + quadrupole 
(squares). The results are shown for modes of degrees I = 1 (left) and I = 2 (right). In the case of a magnetic field composed of a 
dipolar plus a dipolar component both components have the same polar strength, which is such as to make the total polar strength of 
the magnetic field is B p =3 kG. 
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Figure 5. Real part of the phase shifts as function of co-latitude 
for two modes of consecutive radial order (and degree 1 = 1) 
in a model permeated by a quadrupolar magnetic field of polar 
strength b q = 3kG. The upper panel shows the results for the 
modes with frequencies immediately below (continuous line) and 
immediately above (dashed line) the first frequency jump while 
the lower panel shows the results for modes with frequencies im- 
mediately below (continuous line) and immediately above (dashed 
line) the second frequency jump. 
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Figure 6. Same as figure[S]but for a magnetic field which configu- 
ration is the sum of aligned dipolar and quadrupolar components 
with the same polar strengths. 



a cyclic behaviour with the magnetic fiel d intensity should 
also b e expected. In their stellar models ISaio fc Gautschvl 
( 2004) found both a scaling with magnetic field intensity and 
oscillation frequency, and a cyclic behaviour of the frequency 
shifts. The scaling the authors found for the frequency shifts 
in their models, namely, f{B^' 7 v), differs, in the exponent, 
from the scaling found in polytropes. As explained in their 
paper, this difference reflects the difference in the structure 
of the outer layers of the two models. 

In the present models the frequency shifts also display 
an approximate cyclic behaviour, which may be seen in fig- 
ure Q A scaling of the frequency perturbations with mag- 
netic field and oscillation frequency, of the type f{Bp A v), is 
also evident in the same figure. The exponent of the scaling 



is different fro m both the exponent fo und in the polytropic 
models and in lSaio fc Gautschvl (I2004I) . which was to be ex- 
pected given that the outer layers of the model used here 
are different from the outer layers of both those models. 

We should note that the scaling shown in figure [7| seems 
to be valid only in the range of magnetic field intensities 
shown, namely, 3-5 kG. To adjust the frequency shifts ob- 
tained with less intense magnetic fields to the curve shown 
in figurc|7|we would need to change the exponent in the scal- 
in g; again. This is not surp rising, however, since as argued 
bv lSaio fc Gautschvl j2004l) the exponent in the scaling re- 
flects the average structure of the magnetic boundary layer. 
The extent of the magnetic boundary layer depends on the 
magnetic field strength (reaching deeper layers if the mag- 
netic field is more intense). Thus, it is not surprising that the 
scaling with the magnetic field (with fixed exponent) breaks 
down when the range of magnetic field intensities considered 
gets too wide. 



5.4 Comparison with Saio's results 

To compare our results with those obtained with Saio's code 
it is important to eliminate the effects brought about by the 
differences in t he models used. Since the outer layers of the 
model used by ISaio fc Gautschvl i2004l) were significantly 
different from those of the model used in this work, Saio 
(private communication) has computed new frequency shifts 
using the model described in tableQwith a dipolar magnetic 
field. A comparison between the results obtained with the 
two codes is shown in figure [S] 

Although the results computed with Saio's code for the 
model used here are qualitatively in agreement with our re- 
sults, there are significant differences in the details. On the 
lowest left hand side panel we can see that the dependence 
of the imaginary part of the frequency shifts on the degree of 
the modes is consistent for both codes (with I = 1 showing 
the largest shifts and I = 2 the smallest). However, the same 
agreement does not hold for the real part of the frequency 
shifts. On the right hand side panels we can also see that the 
frequency shifts computed with the two codes jump at differ- 
ent frequencies. This might result from the differences in the 
boundary condition applied at the top in the two codes. But 
probably the most striking difference between the two sets of 
results is the amount by which the real part of the frequency 
shifts jump, which is significantly l arger in our result s . This 
difference was already stressed in ISaio fc Gautschvl j2004l) 
for the case of polytropic models. Curiously enough, the fre- 
quency shifts seem to come into agreement again when the 
unperturbed frequency is sufficiently far from the jump. 



6 CONCLUSIONS 

In the present paper we investigate on the effect that mag- 
netic fields with different configurations have on the oscil- 
lation frequencies of acoustic waves propagating in the in- 
terior of roAp stars. We consider both magnetic field con- 
figurations which produce perturbations that are symmet- 
ric about the magnetic equator (like pure dipolar and pure 
qudrupolar magnetic fields) and a 'mixed' magnetic field 
configuration, composed by the sum of aligned dipolar and 
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Figure 7. Scaling of the real (top) and imaginary (bottom) parts of the frequency shifts with magnetic strength, according to the relation 
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symbols) and Cunha's code (filled symbols), for the stellar model described in table HI and the same magnetic field, with dipolar 
configuration and polar strength of B p = lkG. The left panels show the comparison for modes of different degrees I = (circles), 
I = 1 (triangles), 1 = 2 (squares), 1 = 3 (diamonds) before the first frequency jump. The right panels show the compassion for modes of 
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quadrupolar components, which induces a magnetic pertur- 
bation that is not symmetric about the magnetic equator. 
The interest of studying magnetic fields that may produce 
perturbations that are not symmetric about the magnetic 
equator may be exemplified with the well known roAp star 
HR 1217. The modes observed in this star are believed to be 
of alternating even and odd degree. However, when attempt- 
ing to interpret the multiplet structures seen in the power 
spectrum of this star one finds that odd degree components 
need to be included when describing any of the modes as 
a sum of spherical harmonics. This distortion of the even 
degree modes into a sum of spherical harmonics that in- 
clude odd components may not be produced by dipolar or 
by quadrupolar magnetic fields. However, it could be pro- 
duced by a magnetic field which effect on pulsations is not 
symmetric about the magnetic equator. 

The main conclusion of the present paper is that the 
magnetic frequency perturbations are influenced in two dis- 
tinct ways, associated with two classes of effects: 



• the overall pattern followed by the frequency shifts 
scales with frequency in a way that depends essentially on 
the structure of the outer layers of the model and on the 
intensity of the magnetic field. 

• the amount by which the real part of the frequency 
shifts jump at well defined frequencies depends essentially 
on the geometry of the problem, i.e. on the magnetic field 
configuration and on the degree of the mode. 



We also confirmed that the latter result is a direct conse- 
quence of the behaviour of the phase shifts calculated at the 
bottom of the magnetic boundary layer, which jump by n 
at well defined values of co-latitude, namely the co-latitudes 
at which \B X \ = 1/2|S Z |. This simple dependence of the 
amount by which the real part of the frequency shifts jump 
on the geometry of the problem opens an interesting possi- 
bility, namely, that of using anomalies that might be found in 
the power spectra of roAp stars to infer information about 
the magnetic field configuration and/or the degree of the 
modes observed. 

In the present work we have also investigated on the 
effect of changing one of the boundary conditions applied 
at the outermost layer of the stellar model. In particular 
we have considered the case in which the acoustic compo- 
nent of the wave that may be identified in the magnetically 
dominated region is allowed to propagate away when its fre- 
quency is above the appropriate acoustic cutoff frequency, 
and the case in which that component of the wave is is fully 
reflected. 

We found that the boundary condition applied to the 
acoustic component of the wave in the outer layers of the star 
influences in a significant way the acoustic wave in the inte- 
rior only when the oscillation frequency is above the acous- 
tic cut-off frequency of an otherwise similar model with no 
magnetic field. Moreover, we found that the differences in 
the results obtained with the two boundary conditions are a 
consequence of the effect the boundary condition has on the 
solution when the magnetic field is only moderately inclined 
in relation to the local vertical coordinate. 



7 ACKNOWLEDGMENTS 

I am very grateful to Hideyuki Saio for providing part of 
the data used in figure |H] MC is supported by FCT (Portu- 
gal) through the grant BPD/8338/2002. This work was sup- 
ported by FCT and FEDER (through POCI2010) through 
theproject POCTI/CTE-AST/57610/2004" 



REFERENCES 

Balmforth N. J., Cunha M. S., Dolez N., Cough D. O., 

Vauclair, S., 2001, MNRAS, 323, 362 
Bigot, L., Provost, J., Berthomieu, C, Dziembowski, W. 

A. & Goode, P. R., 2000, A&A, 356, 218 
Biront, D., Goossens, M., Cousens, A. & Mestel, L., 1982, 

MNRAS, 201, 619 
Campbell, C. & Papaloizou, J., 1986, MNRAS, 220, 577 
Cunha, M.S. ,2001, MNRAS 325, 373 
Cunha, M.S.,2002, MNRAS 333,47 

Cunha, M.S., Fernandes, J.M.M.B. & Monteiro, 

M.J.P.F.G, 2003, MNRAS 343, 831 
Cunha, M.S. ,2003, in M.J.Thompson and J.Christensen- 

Dalsgaard (eds), Stellar Astrophysical Fluid Dynamics. 

Cambridge university Press, UK, p. 51 
Cunha, M.S. ,2005, JAA, in press 

Cunha, M.S. & Cough, D.O., 2000, MNRAS 319, 1020 
Dziembowski, W. A. & Goode, P. R., 1996, ApJ 458, 33 
Elkin, V.G, Riley, J.D., Cunha, M.S., Kurtz, D.W., 

Mathys, G, 2005, MNRAS 358, 665 
Cough, D.O., 1993, Astrophysical fluid dynamics, J-P. Zahn 

& J. Zinn- Justin (eds), Elsevier, Amsterdam, p. 399 
Cough, D.O., 2005, in The Roger Tayler Memorial lectures, 

Astronomy & Geophysics, special issue, Sue Boweler (ed), 

Royal Astronomucal Society (pub.), p. 16 
Hubrig, S., Szeifert, T., Scholler, M., Mathys, G, Kurtz, 

D.W., 2004, A&A 415, 685 
Hubrig, S., Nesvacil, N., Scholler, North, P., Mathys, G, 

Kurtz, D.W., Wolff, B., Szeifert, T., Cunha, M.S., Elkin, 

V.G, 2005, A&A, submitted. 
Kurtz D. W., 1982, MNRAS, 200, 807 
Kurtz D. W., Martinez, P., 2000, Bait A 9, 253 
Kurtz, D. W., Elkin, V.G, Mathys, G, Riley, J., Cunha, 

M.S., Shibahashi, H., Kambe, E., 2000, BaltA 9, 253 
Leckrone, D.S., 1973, ApJ, 185, 577 
Mathys G. & Hubrig S., 1997, A&AS, 124, 475 
Matthews, J.M., Kurtz, D.W., Martinez, P., 1999, ApJ 511, 

422 

Morel, P., 1997, A&AS 124, 597 

Priest, E.R., 1982, Solar Magnetohydrodynamics, Dor- 
drecht: D. Reidel Publishing Company 
Roberts, P. H. & Soward, A., 1983, MNRAS, 205, 1171 
Saio, H. & Gautschy, A., 2004, MNRAS 350, 485 
Saio, H., 2005, MNRAS 360, 1022 

Stepien, K., 1994, in J. Zverko, and J. Ziznovsky (eds.), 
Chemically peculiar and magnetic stars, Astronomical In- 
stitute, Slovak Academy of Sciences, pg. 8 



Pulsations in magnetic Ap stars 11 



APPENDIX A: BOUNDARY-LAYER EQUATIONS 

In the magnetic boundary layer we adopt the approximation of a locally uniform (inclined) field, which is essentially force 
free and which therefore does not influence the stratification of the surface layer. Thus, from equations Q-JSJ , we find that 
in that layer low-amplitude adiabatic magnetically nondiffusive pulsations are described by the following system of linearized 
equations: 

J| = -Iv Pl + 4Vp+— (VABi)AB, (Al) 
ot p pop 

pi = -€ • Vp - 7PV ■ £, (A2) 
Pi = -t ■ Vp - pV ■ €, (A3) 

Bi=VA(£AB), (A4) 

where the subscript 1 denotes an Eulerian perturbation and p, p, and B, now stand for the corresponding unperturbed 
quantities. 

We adopt a plane-parallel approximation at each latitude and longitude of the star. The local vertical coordinate, z, is 
taken to increase outwards, and to be zero at the surface of the star. The horizontal coordinates, x and y, are chosen to form a 
right-handed locally Cartesian system with the z axis, orientated such that the x axis is parallel to the horizontal component 
of the magnetic field. Thus we set B = (B x , 0, B z ). Moreover, in accordance with our assumption of a locally uniform magnetic 
field, we neglect the derivatives of B x and B z in equations IAH — (IA4t . 

We seek solutions of the system of equations HA1I - IA4I that depend on x, y and time t through the factor 
exp [i (k x x + k y y + u)t)], where k= (k x ,k y ,Q) is the (local) horizontal wave number and u> is the oscillation frequency, k 
can be related to the degree I and the azimut hal order m of the mode co nsidered by means of a standard asymptotic ap- 
proximation for the spherical harmonics fe.g.. lCampbell fc Papaloizoulll98^1 . For axisymmetric modes the relation reduces to 
k x ~ L 2 /R and k y = 0, where L = I + 1/2. 

We decompose the displacement into its vertical component and its horizontal components parallel and perpendicular to 
k, defined respectively by, £ z = £.e z , u = £.k/|k| and v = £. (e z A k) /|k|, where e z is a unit vector in the z direction. Using 
equations JA2I — iA4t to eliminate pi,pi and Bi from equation 1A1L and neglecting the derivatives of B x and B z , we obtain 

- tu 2 pu = i\k\W + — (B ■ Vfu - ^^B ■ V (V ■ £) , (A5) 

Mo A*o|k| 

= — (B-V)W^fB-V(V-O, (A6) 
Mo /Uo|k| 

- Lo 2 pi z = ^- 9 V- [pQ — — [B • V (V • £)] + -(B • V) 2 ^ z , (A7) 

az p,o p,o 

where, 

W = i-Vp+ (~/p+— ) V-£-— B-V(B-g) (A8) 
V Mo / Mo 

and 

(A9) 

p az 

Next we express the boundary-layer equations in a dimensionless form, by defining 

Z U3 p p U 

V — n> O" = — , P=— , Q=—, U =S> 

H cro po po H 

v ^ z Bi 

V — — , Z = — , K x — Rk x , K y = Rk y , bi = - - 1/2 , 

H H (popo&o 2 R 2 ) ' 

where 7? is the stellar radius, ao = (GM / R 3 ) 1 ^ 2 , po and po are, respectively, the central pressure and density and u, v and £ z 
are the z-depend parts of the corresponding components of the displacement defined by ^u, fi,^ z j e 1 { k ^ x + k vy+ bjt ) — (u^^z). 
Introducing the new variables in the system of equations (IA5II — (IA9I) . and denning 
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we obtain 

b x b g ^Z" + i \CypK + bl?f \ Z' + (iCp'K - b x b z KK x ) Z + ib x b z K y V' + b 2 x K x K y V 
-b 2 z U" + (-a 2 g + CjpK 2 + b 2 K 2 - b 2 x K 2 ) U = 0, 

. - (C 7P + b 2 x ) Z" - Cyp'Z' + [-o 2 e + b 2 x K 2 x +C - p")) Z - b x b z ^V" + ibl^V + b x b z ^U" 
+i UypK + bl?f \ U' + [iC ((7 - 1) p'K - j'pK) - b x b z KK x ] [7 = 0, 

. -b x b z ^Z" + ibl^^Z' - b 2 z V" + 2ib x b z K x V + (-a 2 g + b 2 x K 2 ) V + ib x b z K y U' + b 2 x K x K y U = 0, 

where prime denotes differentiation with respect to r\ and K — (K 2 + K 2 ^ 1 ^ 2 . For axisymmetric modes, considered in this 
work, K y — and K x = K. Moreover, in this case the component V of the dimensionless displacement decouples from the 
components Z and U , and the third-order system of second-order complex differential equations given above separates into a 
second-order system describing the components Z and U and a single second-order equation describing the V component. 

The second-order system of second-order complex differential equations describing the components U and V was then 
transformed into a fourth-order system of first-order complex differential equations, and integrated using a fourth-order- 
accuracy Runge-Kutta algorithm. 



APPENDIX B: BOUNDARY CONDITIONS 
Surface boundary conditions 

Sufficiently high in the atmosphere, where the magnetic pressure is much larger than the gas pressure, the response to motions 
that are perpendicular to the magnetic field is dominated by the perturbed Lorentz force. Motions along the magnetic field 
lines, however, will generate an acoustic response, associated with the gradient of the perturbed gas pressure. Hence, in this 
region two decoupled waves may be identified: a compressional Alfven wave, which displacement is essentially perpendicular 
to the direction of the the magnetic field, and an acoustic wave, which displacement is essentially along the latter. 

Because the density tends to zero as we move towards the outer layers of the atmosphere, the magnetic field must tend 
to a vacuum field implying that: 

(VABi), = 0, (Bl) 
or equivalently, 

[(B ■ V){ - BV ■ €]. = (VVm)s, (B2) 
where tpm is an external magnetic potential, regular at infinity and satisfying Laplace's equation, 

V>m = Ae'^ z e i ^ x+kyV+ut \ (B3) 

A is a constant and the subscript s indicates that the condition is to be applied at the surface. 

The matching into a vacuum field is translated into a condition on the compressional Alfven wave. Combining equations 
l|B2fl and l|B3|l . and writing the result in terms of dimensionless variables, we find, for axisymmetric modes, 

(e ± ' - Ke ± ) B = 0, (B4) 

where e± = Z cos a — U sin a is the component of the dimensionless displacement perpendicular to the direction of the 
magnetic field and a is the angle between the magnetic field direction and the x axis. This boundary condition is applied 
sufficiently high in the atmosphere, in a region where e' L / e x 1/H (which implies that (aH/b) 2 g <C 1), with H = p'/p. 

Next we consider two possibilities for the acoustic wave associated with the displacement parallel to the direction of the 
magnetic field. In some models this component of the solution is artificially reflected by imposing the condition 

(V • Os = 0, (B5) 

which, in terms of dimensionless variables becomes, 

(Z' - iKU) B = 0. (B6) 

In other models the same component of the solution was matched into the approximate analytic solution appropriate for an 
isothermal atmosphere. Starting from equations (IA5ft — (IA9II adapted for axisymmetric modes in an isothermal atmosphere, 
and considering the limit when the magnetic pressure is much larger than the gas pressure, we find that the dimensionless 
component of the solution parallel to the magnetic field direction, ey = Z sin a + U cos a, obeys approximately the equation 

"^f 1 9- K \ ' 1 ( r2 ■ g ^ n rtvf\ 

e ll + 77 - 2i— — U|| + TT—njfZ - 7™n: - t -rjr— U|| = 0. (B7) 



H tan ay \C7psh1 a tan a _fftana 
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The former equation accepts the solution 

e„ = A e K^, + « (B8 ) 
where A is a constant, <j> is a phase and 



Ka ~ tan a ± Y C 7P sin 2 a 4ff 2 ' (B9) 

The numerical solution corresponding to the displacement parallel to the magnetic field direction is matched at the surface 
onto the analytical solution 1B8L by imposing the condition, 

~ + iK a ]e {] =o. (BIO) 

When the argument of the square root is positive, K a is real and the sign in equation 1B91 is chosen such as to make the 
solution correspond to an outwardly propagating wave. When the argument of the square root is negative, K a is complex and 
the sign in equation 1B9I is chosen such as to make the solution correspond to an evanescent wave. 



Interior boundary condition 



In the interior, where the magnetic pressure becomes negligible when compared with the gas pressure, the magnetoacoustic 
wave decouples again into a wave which is essential ly magnetic, a slow Alfv en wave, and a wave which is essentially acoustic. 
Here we follow the boundary condition applied bv lCunha fc Gouehl i2000T) . To leading order the acoustic component in the 
interior of the star is described by the system of equations HA5II - 1A9II with B set equal to zero, namely, 



u pu a = — i|k|pi, 



and 



2 , dpi 
u p£ za = — +gpi- 
oz 



(Bll) 



(B12) 



Moreover, the displacement associated with the slow Alfvenic compon ent is essentially along the horizontal direction (along 
the x axis for axisymmetric modes). Following lRoberts fc SowardHl983l) . this component is assumed to dissipate in the interior 
and thus the corresponding numerical solution is matched onto the analytical asymptotic solution for an inwardly propagating 
slow Alfven wave, derived in the JWKB approximation, namely, 



Dp 



-1/4 



exp 



\ B 2 



1/2 



dz 



.kB x 



exp [i (kx + cot)] , 



(B13) 



where D is a (complex) constant and u m = u — (i|k|pi)/(w 2 p) is the magnetic part of the horizontal displacement. In terms 
of dimensionless variables, the matching condition thus become, 



with 

t/ m = 



a 2 g 
7,2 

17- 



1/2 



Oz 



u 



.CK 



C-ypK 2 



(Zp' + jpZ') 



(B14) 



(B15) 



The boundary conditions defined above allow for the complete determination of the solution of the equations in the 
magnetic boundary layer (except for the amplitude which cannot be determined under the linear approximation). Thus, 
we may then use the numerical solution for the vertical component of the acoustic mode in the interior to match onto the 
corresponding asymptotic solution. Using equation 0, that matching is translated to the condition (in dimensionless variables 
defined in the plane-parallel approximation), 



z 



— ~ — /Ctan<5, 



(B16) 



which is applied at the matching depth tj* = 1 — R* / R, with K, = Rk. Equation <B16> allows us to calculate the complex 
phase 8 up to a factor of n. This paper has been typeset from a TeX/ I^TjgX file prepared by the author. 



